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Abstract 



Using the supersymmetry method we analytically calculate the local den- 
sity of states, the localiztion length, the generalized inverse participation ra- 
tios, and the distribution function of eigenvector components for the superpo- 
sition of a random band matrix with a strongly fluctuating diagonal matrix. 
In this way we extend previously known results for ordinary band matrices 
to the class of random band matrices with preferential basis. Our analytical 
results are in good agreement with (but more general than) recent numerical 
findings by Jacquod and Shepelyansky. 
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Typeset using REVTgX 



Banded matrices with random elements play an important role in the description of both 
certain systems exhibiting quantum chaos (in particular the kicked rotator 0) and electron 
transport in quasi Id geometries ||. In a series of papers [|3],|J, Fyodorov and Mirlin have 
performed a detailed analytical investigation of the properties of such matrices, making use 
of Efetov's supersymmetry technique || . In recent times particular interest in random band 
matrices (RBM) with preferential basis (PB) (realized, e.g., by a very strongly fluctuating 
diagonal) has emerged for at least two different reasons: First, any attempt to go beyond 
the quasi Id case for the electron transport problem, either by employing a Fokker-Planck 
approach |J or the (equivalent [|7|) a model formulation necessarily leads away from 

the isotropic situation and introduces a model with preferential basis. Second, Shepelyansky 
|| recently studied the problem of two interacting particles in a Id random potential by 
reducing the Hamiltonian to a RBM with PB. He found that the two-particle localization 
length is strongly enhanced as compared to the one-particle localization length. This result 
was reinforced, made more precise and further investigated by subsequent work (T^jni • Very 
recently, Jacquod and Shepelyansky |L2] studied numerically certain properties of RBM with 
PB, especially the local density of states, the inverse participation ratio and the level spacing 
statistics. 

In this letter, we extend the analytical treatment of Fyodorov and Mirlin |3],|]] to the case 
of RBM with PB, derive explicit formulas for the local density of states and the localization 
length, and present expressions for the generalized inverse participation ratios and the dis- 
tribution function of eigenvector components in terms of previously known results for RBM. 



Our results are in good agreement with the asymptotic estimates given in |L2 and can also 



explain certain numerical deviations which occur in [0 in the limited range of accessible 
system parameters. 

We consider the random matrix 

Hij = rjijSij + Cij, (hj = !, • • • , N ), (!) 

where the 7/j = rju are real random numbers with the distribution function po{v)- We 
introduce a scale parameter Wb by setting (jf) « W b 2 . The matrix ( is either symmetric 
{Cij = Cji-> P = 1) or Hermitian (Qj = Q { , (3 = 2) with Gaussian random variables satisfying 
(\0j P) = (1 + <%(2 — /3))Aj/2. Aij — a(\i — j\) is a function that decays on the scale of the 
bandwidth b and has typical values l/y/b. In other words, ( is exactly the RBM considered 
earlier by Fyodorov and Mirlin [H] and rj introduces a PB. For later use we define (as in M) 
Bo = J2r a ( r ) an d B 2 = X)r a ( r ) r2 /2- The class of matrices introduced in (P contains the 
case considered in as a particular example: 

a(r) = 57TTl5 e(6 - |r|) - (2) 

The importance of the r\i is governed by the ratio Wb/\/b of the spacing Ar/ ps Wb/b of those 
r/i which are coupled and the typical value Ct yP . ~ b~ 1 / 2 of the coupling matrix elements. We 
can distinguish three important regimes: (i) Wb/Vb 3> 1 =>• the coupling matrix elements are 
very weak and can be treated perturbatively, (ii) 1/Vb Wb/ Vb <C 1 =>- the r/i are still much 
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larger than the Qj : but the coupling to the Qj becomes nontrivial, and (iii) Wb/Vb ^$1/ \fb =r- 
the rji are comparable to the Qj and the RBM results in |^,f|] are applicable. 

In the regime (i), one may apply simple perturbation theory to calculate the eigenfunc- 
tions. This yields the behavior \ip(jb)\ ~ |^(0)| (Ctyp./Ar/)-?. The corresponding localization 
length is then estimated as £ ~ 2b/ \n(W b 2 /b), a behavior that was also anticipated (and 
numerically confirmed) in || . In the present work, we are mostly concerned with regime (ii) 
and the crossover to (iii). 

We are interested in calculating the (position dependent) generalized inverse participation 
ratios (we use the notation of ||): 

P,(E, n) = IT Wn)\ 2 'S(E - &)) 

= !-- m , gfo" ' " (7+ 't 2)?" C^»'<^'"') <«= i+m >- < 3 > 

To perform the average over the product of Green's functions we use the supersymmetry 
method. For all details and the standard notation we have to refer the reader to P,[T3"|]. 
We define a supersymmetric functional 

F(J) = (sdet~ 1/2 (£ - H + isA + J)) = (sdet~ 1/2 (l + gj) 



l + x+G+J-^l + x-G^,)- 1 ^, (4) 

where g = diag(G + , G~), and J = diag(x+Pe, x_Pb) ® e^e^ = J (g> e n e|,. Here, x + and x_ 
are simple numbers, Pb is a projector on bosonic variables, e n is a vector with (e n )i = 5 n i, 
and (. . .)( denotes averaging over (. Our formulas are valid for both orthogonal and unitary 
symmetry, provided the four-dimensional graded space in the unitary case is simply doubled. 
Following standard procedures Hf|], we express @ as an integral over supervectors, average 
over Qj, perform a Hubbard-Stratonovitch transformation and finally arrive at 



F(J) = J L>[<7i]exp{-A(<T)}, 



= o J2 str(a i (A- 1 ) ii a j ) + str ln ( E ~ Vj + ieA + °il 2 + J6 Jn), (5) 
8 ij 2 i 

where the <7j are 8x8 supermatrices. Note that the functional still depends on the particular 
realization of the rjj. No average over these variables has been performed yet. The saddle- 
point condition for (^) reads 

<71 = -Y A H^ —* T- ■ ( 6 ) 

With the ansatz a = To + iY\Q for the solution of (0), where Q satisfies Q 2 = 1 and is a 
proper element of the coset space introduced in ||, we have for z = T + iTi. 

z = -B J d VP o(v) E _^ + z/2 (M*0 > 0). (7) 



3 



Here we have replaced each element of the sum in (^|) by its average over 77. This procedure 
is only justified in the limit of "overlapping resonances" , i.e. if the widths 1^ of the super- 
imposed Lorentzians are larger than their spacings determined by those T]i contributing to 
the sum. With the result Ti ~ p(E) ~ 1/Wb (see @ below) this immediately leads back to 
the condition Wb <C V~b, which is fulfilled in our regimes (ii) and (iii) defined above. Also, in 
precisely this limit we can expect the saddle-point solving @ to be homogeneous in space 
as in our ansatz. 

Equation (^) determines T and r x and hence the local density of states (p(j, E))^. Using 
(|) it is not difficult to show that (G^)^ = (E — r)j + z/2)^ 1 and therefore 

w, E )}( = -\i m{ a mk = (g + r ^ )a + n/4 . (8) 



This is exactly the Breit-Wigner form found numerically in ||12|| . Averaging this expression 
over all sites (assuming that N^ 1 ■■) = (■■ -)r?) and using the saddle-point equation (|7|) 

we get 

p(E) = ((p(j,E)) c ) v = V 1 /7rB (9) 

for the average density of states. From (0) and (^), we find the limiting cases p(E) ~ po(E) 
for W b > 1 and p(E) ~ y/2B - E" 2 /(TtB ) for W b < 1. 

Inserting the saddle-point solution in (K|) we proceed in analogy to MM and arrive at 



F(J) = J D[Q i ]exp{-£ 2 (Q)}, 



C 2 (Q) = E ^[(Qj+t - Qjf] + \ E str M 1 + MjiQ^ieA + J5 jn )l 
M d (Q) — (E — m + T /2 + iT 1 /2Q)-\ (10) 
Here, £ = 2Y\B 2 /Bl = 27r 2 p(E) 2 B 2 denotes the localization length |14]J5|1 for the wave 



functions fI5|| , |-0(ri) | ~ \ip(0)\ exp(— n/£). Due to the one-dimensional structure of the 
Q-functional (0) F(J) can be written as 

F{J) = JdQ Y(Q, n) Y(Q, N - n) sdef-^l + M n {Q)J), (11) 

where the function Y(Q,n) arises from integrating out all Q-matrices up to site n. The 
properties of these functions have been discussed in detail in BM. From (f|) it is clear that 

\l+m 

((GL) l (G- n ) m ) = L_i_^5m jr(j)| B±=o . (12) 

Performing the source term derivatives in (|11|) and using that 

_ E — r/j + r /2 - zgV2)Q 

AQ) ~ (e — r]j + r ) 2 + rf/4 (13) 

we finally get in analogy to [|J (for e — > and (3 = 2) 



4 



((G+J(G- n ) m ) = [ + m j {-m P C{ Vn )) l+m J dQY(Q,n)Y(Q,N -n)Q[ liBB Q^ BB (14) 
with 

1 Ti/2 
C{7]n) = ~p(E- Vn + T /2) 2 + T 2 /A (15) 

Comparing (|T4]) with the corresponding expressions in 0] we can immediately conclude that 

P q (E,n) = C(r ]n yP q FM (E,n). (16) 

The superscript "FM" refers to the result by Mirlin and Fyodorov Q for ordinary RBM, 
where their rescaled length x has to be identified with x = 2N/(/3£) [£ as in eq. fllO)]. 
Similarly, following the steps in Q to derive the local distribution function V(y, n) for y = 
iV|^(n)| 2 defined by P q (E,n) = N 1 ~ q J °° dy V(y, n) y q , we find M 



rM = cK) pFM {cbY n )- (17) 

We see that the local distribution is rescaled by the energy dependent factor C(rj n ). In the 
resonant case, E ~ r] n , the typical amplitude y exhibits a strong enhancement (in the limit 
Wb S> 1, i.e. Ti ~ W^ 1 ), whereas for (E — r] n ) 2 ^> B a strong suppression is observed. 

Up to now we have considered a fixed realization of the diagonal elements rjj. The 
additional average over the sites (or equivalently over 7] n ) results in 

p « = h»u> ck v "'{ck)- (l8) 

This completes our formal derivations for the case of RBM with PB. In the following we 
derive some more explicit results for the special case of Shepelyansky (see (0)). 

With B = 2/3 we have p(E) = 3T 1 (E)/2ir, where T^E) is determined by (0) and 
depends on E and po(rj) (i-e. on Wb for the case (|2|)). In the main part of Fig.l we show 
p(E) for different values of Wb, demonstrating the crossover from a semicircle to the box 
shape (0). For E = we must have r = and (|7|) simplifies to 

T l (0) = — arctan [ -^-) . (19) 
u ; SW b V r i(°)/ 

With E> 2 = b 2 /9 we have £ = r!(0) 2 5 2 /2 so that we can express £ as a function of Wb- We 
recall that our treatment is valid for all Wb with < Wb <C Vb, i.e. for RBM, RBM with 
PB and the whole crossover. In the inset of Fig.l we have plotted £ versus Wb- Interestingly, 
the crossover region between the two limiting cases £/b 2 = const. (RBM) and £/& 2 ~ ^/W b 2 
(RBM with PB as specified in (0)) is a comparatively small interval around Wb = 1. 

For the asymptotic case Wb ^> 1 considered in the estimates in P,|T2"[ we get by expanding 

(0) 
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This is in very good agreement with the estimate b 2 /2W 2 in PJjjJ- We also remark that a 
naive application of the RBM results in || (where one simply modifies the function a(\r\) 
to account for the large diagonals) does reproduce the parameter dependence in (pD|) but 
yields a (wrong) prefactor 4/3 instead of 7r 2 /18. 

Finally, we calculate the quantity £,ipr, defined by £Jpr — P% — N^ 1 J2 n P2(E, n) in the 
limit W b 1. We get for arbitrary £ and N 

vr 2 / 3 lr 1 
e/Pi? -T2W \2N + e ' (21) 



with the metallic limit (N -C £) 



n 2 N N 

£ IPR = « 0.548^, (22) 

18 W^ 2 



and the localized limit (N £) 



«"« = ilw 5 = £1? M 0451 4 (23) 



These analytical results compare quite favourably with the estimates N/2W 2 (resp. b 2 /AW b ) 
in [0. Also, in view of the exact formula (|2T|), which interpolates between the metallic and 
the localized regime, the numerical deviations from £,ipr ~ W b 2 (resp. W b A ) found in [12| 
for non-asymptotic system parameters come as no surprise. 

Concerning the distribution (|18|), let us briefly discuss the simplest case (with /3 — 2 and 
£ ^> AT), where V FM (y) = exp(-y) as for a GUE random matrix. At E — and W& ^> 1 
[Ti = 7r/(3W 6 )], we find: 

1 r w b 

V ^ = W / ^3(^ + r?/4) exp[-3(^ + r?/4)y]. (24) 

The probability is shifted to very small amplitudes with y <C W b ~ 2 , i.e. V(y) ~ W 2 , and 
also to higher amplitudes y 3> W b ~ 2 , V(y) ~ exp[— 7r 2 ?//(36W 6 2 )]. 

In conclusion we have applied the supersymmetry method to the class of RBM with PB. 
The key to technical progress was the derivation of a a model for a fixed realization of the 
strongly fluctuating diagonal matrix elements. The average over these diagonal variables 
was then performed in a second step, whenever appropriate. Our results account for the 
complete crossover between ordinary RBM (as treated earlier in |||4[]) and RBM with strong 
PB. For the specific model chosen in recent numerical calculations |12] the agreement with 
our findings is very satisfactory. In particular, the existence of two different scales charac- 
terizing the wavefunction as evidenced by the difference between £/p# and £ in fl2~l"|) has been 
demonstrated analytically. We believe that further investigations of RBM with PB along 
the lines given in this letter are interesting and feasible. 

We are grateful to J.-L. Pichard and D. Shepelyansky for instructive discussions. This 
work was supported by the European HCM program (KF) and a NATO fellowship through 
the DAAD (AMG). Upon completion of this manuscript a preprint with very similar results 
by Y.V.Fyodorov and A.D.Mirlin ( |cond-mat/9507043|) came to our attention. 
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FIG. 1. The average density of states for the particular band matrix considered in Refs. [|,12] 
(see Eq. (|2|)) and the values Wb = 0, 1, 2, 10. The inset shows the localization length £ normalized 
by b 2 at E = as a function of Wb (in a doubly logarithmic representation). 



S 



